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Abstract 

We study the photo-association of Bose-Einstein condensed atoms into molecules using an op- 
tical cavity field. The driven cavity field introduces a new dynamical degree of freedom into the 
photoassociation process, whose role in determining the stationary behavior has not previously 
been considered. The semiclassical stationary solutions for the atom and molecules as well as the 
intracavity field are found and their stability and scaling properties are determined in terms of ex- 
perimentally controllable parameters including driving amplitude of the cavity and the nonlinear 
interactions between atoms and molecules. For weak cavity driving, we find a bifurcation in the 
atom and molecule number occurs that signals a transition from a stable steady state to nonlinear 
Rabi oscillations. For a strongly driven cavity, there exists bistability in the atom and molecule 
number. 
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I. INTRODUCTION 



Molecules can not be directly cooled using the laser cooling techniques that led to Bose- 
Einstein condensation of alkali atoms because of the complex rotational- vibrational spectrum 
of the molecules. As a result, two-color Raman photoassociation and Feshbach resonances 
have emerged as tools to create translationally cold molecules starting from ultra-cold atomic 
gases. The conversion of a macroscopic number of quantum degenerate atoms into molecular 
dimers starting from either a Bose-Einstein condensate (BEC) [l[ [3, [3, [4] or a Fermi gas 



8| has been observed by several experimental groups using Feshbach resonances. 



This work culminated in the formation of a molecular Bose-Einstein condensate (MBEC) 

Although Feshbach resonances have been the most successful tool for creating quan- 
tum degenerate gases of molecules, experiments have demonstrated that two-color Ra- 
man photoassociation can also be used to create molecules in the electronic ground state 
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141 ]. Two-color photoassociation has the added benefit that the frequency differ- 



ence between the two optical fields can be used to select a particular rotational-vibrational 
state [15I, Q- This gives photoassociation a potential advantage over Feshbach resonances 
since photoassociation can be used to prepare molecules in their rotational-vibrational 
ground state. While the molecules created via a Feshbach resonance are translationally 
very cold, they are vibrationally very hot and can decay to lower lying vibrational states 
via exoergic inelastic collisions with atoms or other molecules. For an atomic BEC, the 
molecular two-body decay rates are of the order 10~ n — 10~ 10 cm 3 /s, which gives a lifetime 
of 100/is for typical atomic densities 

Hal- 



From the perspective of atom optics, the conversion of atoms into molecules via a Fesh- 
bach resonance or photoassociation is the matter-wave analog of second harmonic generation 
of photons in a nonlinear crystal with a x susceptibility, which has been used to create 
entangled photon states. Recent experiments have shown the phase coherent and momen- 
tum conserving nature of matter- wave second harmonic generation 19|] . However, implicit in 
these analogies to nonlinear optics is that the coupling strength between atoms and molecules 
is not itself a quantum field. 

Here we address the issue of two-photon Raman photoassociation of an atomic BEC 
inside of an optical cavity. In this case one of the optical fields used to induce the atom- 
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molecule conversion is a quantized mode of a driven optical resonator while the other field 
is a laser with sufficient intensity to be treated as a 'classical' undepleted pump. This 
is reminiscent of early work done by Moore and Meystre 20| studying the coupling of a 
zero momentum BEC atoms to atoms with momentum ±K due to a Raman transition 
involving a quantized ring cavity mode and classical pump. Unlike their system, our model 
involves the interaction of four particles: two atoms are 'destroyed' and a molecule and cavity 
photon are 'created' and vice versa. Consequently, the atom-molecule-cavity interaction is 
analogous to x susceptibility in nonlinear optics, which is known to give rise to four- 
wave mixing. Coherent photoassociation inside of cavity therefore offers the prospect of 
novel nonlinear dynamics between the atomic, molecular, and cavity fields as well as the 
possibility to entangle individual photons with completely different chemical compounds- 
atoms and molecular dimers. 

In fact, only a few papers have previously considered photoassociation inside of a cavity 
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22j. However, unlike our model, theirs was based on single photon photoassociation, 



which is impractical for observing coherent atom-molecule dynamics because the molecules 
created are in electronic excited states and can rapidly decay due to spontaneous emission. 



The authors of Ref. 2l|,|22] employed a positive-P distribution originally developed by 



Gardiner and Drummond [23j to analyze the quantum mechanicial 'phase space' dynamics of 
the three coupled bosonic fields. Since their interaction Hamiltonian required the absorption 
of a photon to create a molecule, the equation of motion for the P-distribution involved 
third order derivatives. Such derivatives meant that the resulting equations of motion could 
not be identified as a Fokker-Planck equation, which can be readily solved using standard 



techniques 



231 ] . As we will show here, Raman photoassociation in which a cavity photon 



is created along with the molecule leads to an equation of motion for the P-distribution 
involving only first and second derivatives that is in the form of a Fokker-Planck equation. 
Consequently, intracavity Raman photoassociation leads not only to more stable molecules 
but also equations of motion for quasi-phase space distributions that are more easily solved 
using known techniques. 

Here, we analyze the semiclassical behavior for the atomic, molecular, and cavity mean 
fields. Our goal is to explore the stationary solutions of the coupled nonlinear mean field 
equations as well as the stability conditions for these solutions in terms of experimentally 
controllable parameters. When the cavity is only weakly driven by an external source, the 
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stationary intracavity field remains close to zero while the molecules exhibit a saddle point 
bifurcation in the molecule number (or equivalently, atoms due to number conservation). Be- 
low a critical driving amplitude, there exists both unstable and stable stationary solutions 
while above the critical point there are no stationary solutions but instead only nonlinear 
Rabi oscillations. For strong cavity driving, the system exhibits bistability in the number of 
molecules. The paper is organized as follows: In section II, we present our model for cavity 
assisted photoassociation and derive the equation of motion for the Positive-P representa- 
tion. From the resulting Fokker-Planck equation we obtain equations of motions for the 
expectation values of the fields. In section III, we analyze the solutions of these mean-field 
equations. In section IV, we discuss the implications of our results. 

II. MODEL 

We imagine that we start with a BEC of atoms inside of an optical cavity, as depicted in 
Fig. [TJ The atoms as well as the molecules formed from them can be trapped inside of the 
cavity using a far-off resonant optical trap similar to what has been recently demonstrated 
with single atoms in a cavity [241 ] . At temperatures T « 0, we can assume that all of the 
atoms are in the ground state of the trapping potential with wave function ip a {r). Addition- 
ally, the atoms are assumed to have all been prepared in the same hyperfine state denoted by 
\a). Pairs of atoms in \a) are coupled to electronically excited molecular states |/„), where v 
denotes the vibrational state of the molecule, via a pump laser with Rabi frequency Qi and 
frequency u)\. The pump is treated as a large amplitude undepleted source and therefore 
changes in fli due to absorption or stimulated emission are neglected. 

The excited molecular states are strongly coupled to molecules in their electronic and 
vibrational ground state, \b), via a single cavity mode. Emission of a photon into the cavity 
mode takes a molecule from an excited state to its electronic ground state. Coupling to 
a single mode can be achieved by insuring that only a single cavity mode is close to two- 
photon resonance for the atom-molecule Raman transition and by positioning the atoms and 
molecules around an antinode of the cavity field. The discrete mode structure of the cavity 
allows one to select a particular vibrational state in the electronic ground-state manifold 
of the molecules provided the cavity linewidth, 7, is less than the vibrational level spacing, 
which is on the order of 1GHz [ijj]. This would imply a cavity Q-factor of Q ^> 10 6 , which has 
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FIG. 1: Schematic diagram showing the system under consideration. The kets \a), \b), and \I U ) 
denote pairs of atoms, electronic ground state molecules, and electronically excited molecules, 
respectively, e is the rate at which the cavity is coherently driven by an external laser while 7 is 
the decay rate for photons in the cavity. 



already been achieved with individual atoms trapped inside of a Fabrey-Perot resonator[24]. 
The cavity field frequency is uj c and the vacuum Rabi frequency for the \I V ) — ► \b) transition 
is g cav . Selective emission into a cavity mode offers a distinct advantage over using free 
space spontaneous emission to populate the electronic ground state of the molecules, since 



spontaneous emission would populate a large number of vibrational 



as well as leading to dissociation back into the continuum la, |25|]. This selectivity along 



evels of the ground state 



with the enhanced emission rate in a cavity due to the Purcell effect in the weak coupling 



regime was proposed as a possible matter wave amplifier for molecules 26]. 

The internal energies of states \I V ) and \b) relative to pairs of atoms in \a) are u v and 
/S.U < 0, respectively. We assume that the detuning between the excited states and the 
pump and cavity fields satisfy, 



A v = oj v - ui « (lo u - Au) - u c > \g ca v\, 

where 7" 1 is the lifetime of \I V ) due to spontaneous emission. Under these conditions the 
excited state can be adiabatically eliminated, leading to two-photon Raman transitions 
between \a) and \b) with Rabi frequency 

X(x) w g cav u(x)n*(x)^2ll v I b>u /A u . 



under the assumption that the cavity mode u(x) and Qi(x) are essentially constant over the 
spatial extent of the internuclear wave functions for the molecules and colliding atoms . 
are the Frank-Condon factors for the \£ = a, b) — > \I V ) transitions. 

At zero temperature we can use a single-mode approximation for the molecular field as 
well as for the atoms. The ground state wave function for the center of mass of the molecules 
is denoted by ipb{x)- This approximation for the atomic and molecular field should remain 
valid at least for times long enough to observe dynamics lasting several Rabi periods. For 
longer times, one can expect that some molecules and atoms will be created with excited 
state wave function [28]. The resulting Hamiltonian for the atom-cavity-molecule system is 

Hatom-moi = (Acj - (w, - 0J c ))Pb + Zy (aHe - aW) + h X ad ]2 a 2 + h Xb b ]2 b 2 (1) 

where a, b, and e are bosonic annihilation operators for atoms, ground state molecules, and 
cavity photons, respectively. Moreover, the molecular operators and photon operators have 
been written in a rotating frame, b —>■ bexp[+i(u c — uiijt] and e — > eexp[— iu c t], to remove 
all time dependence from the interaction term. The coupling constant is given by 

ig = J d 3 X1pl(x)x{x)lpa{x). 

We can further simplify things by assuming that the two-photon resonance condition, Au = 
uji — lu c , is satisfied at all times so that the first term in H a tom-moi is identically zero. The 
terms proportional to Xa and Xb represent the two body interactions between pairs of atoms 
and pairs of molecules, respectively. We have not explicitly included two-body interactions 
involving an atom interacting with a molecule because such a term can be written as 



XabN a N b = \xab (n 2 - N 2 - 4N 2 



where iV = N a +2Nb is the total number operator and N a = a) a and A\ = b^b. Since the total 
number of particles is conserved, we are dealing with eigenstates of iV and the atom-molecule 
interaction can be absorbed into a redefinition of the atom-atom and molecule-molecule 
interactions. 

The dynamics of the empty cavity are described by two competing processes. The first 
process is cavity decay, which can be treated using the standard Born-Markov master equa- 



tion for the density operator 



23] 



-^{damping = | (2epe t - e'ep - pe ] e) (2) 



In addition to this, the cavity is pumped by some external source. Coherent driving of 
the cavity by a classical source such as a laser is described by the following interaction 
Hamiltonian, 



pump 



The complete equation of motion for the density operator is then given by, 



dt in P um P ^atom—moli 



p] + — (2epe^ — e^ep — pe*e) 



(3) 



(4) 



The quantum dynamics of the system can be analyzed by deriving an equation of motion 
for the positive-P distribution, P(a, (3, e, a*, (3*, e*, t), which is a representation of the density 
operator in terms of bosonic coherent states, 

\a,(3,e) (a*,/3*,e*\ 



P 



d 2 ad 2 a*d 2 (3d 2 (3*d 2 ed 2 e*P(a, (3, e, a*, f3\ e*, t)- 



(5) 



(a*, P*, e*\a, ft, e) 

where C\a,{3,e) = c\a, (3, e) where C = a,b,e and c = a,{3,e, respectively Note that in 
general a,/3,e and a*,/3*,e* are to be understood as independent complex variables and 
not simply complex conjugates. Using Eq.(jl|) and standard methods 23J], a Fokker-Planck 
equation can be derived for the positive P-distribution P(a, f3,e,t). 





dP 
~dt 



+ 
+ 
+ 
+ 



.da 
d 

da* 

A, 

d(3* y 

d 2 
da* 2 

d 2 



(2ix a O!*a — ga*f3e) 

{-2i Xa aa* 2 - ge*(3*a) + -^-(2i X bi3*f3 2 + ?arV) 



-2ix a aa* 2 + -ea* 2 ) H — 



d_ 

8(3' 
1, d 2 



2t Xa a* 2 + g(3*e*) + 



2 W 
d 2 



-2ixa0? + g(3e) 



dp 2 



2^X ^) /3 2 ) 



2 



-got 2 ) + 



d 2 



i-9<*")) 



P 



(6) 



Un 
in 



ike the equations of motion for single photon cavity assisted photoassociation derived 
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221 ] . there are no third order derivatives with respect to the fields in Eq. [6j Math- 



ematically this originates from the different positions of the cavity field annihilation and 
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22], the interaction was of the 



creation operators in the interaction Hamiltonian. In 
form d 2 b^e + h.c. while here it is d 2 6^ + h.c. 

The most common way to solve Eq. [6] is to map a,(3,e and a*,/3*,e* onto a set of Ito 



stochastic differential equations 



231 ] . Here we are interested only in the dynamics of the 



mean fields, (C) = c, which leads to deterministic c-number differential equations, which 
can be obtained from the Ito stochastic differential equations by setting all noise terms equal 
to zero. The resulting equations of motion are: 

a = —2ix a a*a 2 + ga*/3e (7) 
P = -2i Xb P 2 P* - |« 2 e* (8) 
t = e-\e-\c?P (9) 

The next section analyzes the solutions of Eqs. EHHJ 

III. STATIONARY STATES AND STABILITY 

There exists two distinct cases of stationary solutions in the system. First, there is the 
regime of weak driving, when the stationary value of the cavity field equals zero. The 
opposite case occurs when the driving is strong enough to determine the cavity amplitude 
regardless of the values of the matter fields. 

First we treat the case of weak driving with nonzero stationary values of the atom and 
molecule modes, and rewrite Eqs.([7}|9]) in polar form using 



a(t) = V 'N - 2N b (t)e tea{t) , (10) 
Pit) = y /Ntfje ieb ®- (11) 

Substituting these into Eqs.(J7H9]) we arrive at 

A = A Xa N -2( Xb + ^Xa)N b + g b [ej cos(A) + e R sin(A)] (12) 

2\J l\b 

Nt = -9^fN~b{N - 2N b )[e R cos(A) - e / sin(A)] ) (13) 

e R = e - \e R - \\fW h {N - 2N b ) cos(A) (14) 

ei = ~ej + 9 -^/N b {N - 2N b ) sin(A), (15) 

where e R and ej are the real and imaginary parts of the cavity field mode, and where 
A = Q b — 2Q b . We note that the photo-association dynamics, as given by Eqs.ffT2l)- f|T5l) . 
only depends on the relative phase between atom and molecule modes, A, and the phase 
dynamics of the individual modes is thus of no interest in this situation. 
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As a first case we investigate the possibility of a nontrivial steady state with empty 
cavity e R — ej = 0. From Eq. (fT5l) we see that sin(A) = while Eq. (fT4l) give cos(A) > 0, 
which combined show that A = 0. From Eq. (ll2p . we see that in general, all four equations 
cannot be satisfied if e = 0, but if the nonlinear term is neglected, a steady state is possible. 
Rewriting Eq. (1141) in terms of the molecule fraction, 

z 2 = 2N b /N, 

we find the condition for stationarity 



Pa(z) 



z + C = 0, 



where the dimensionless constant C is given by 



N 



C=-[-r 

9 



-3/2 



(16) 



(17) 



and represents the ratio of the rates for the two competing processes that populate the 
cavity with photons- pumping and emission of a photon during the atom to molecule Raman 
transition. In terms of C, weak driving therefore corresponds to C <C 1 while strong driving 
is C » 1. 

Equation Eq. fflBT) has two roots in the interval < z < 1 if 

2 



c <a 



crit 



3V3' 



(18) 



where the upper limit C cr u denotes the value where min P^(z) = 0, i.e. when the two roots 
coincide, which occurs for z = l/v3. In order to examine the stability of the stationary 
solutions, we linearize Eqs. (jT2l) - (fl"5l) for small perturbations around the stationary values 

A(t) w A + SA(t), (19) 

N b (t)^N + 5N b (t), (20) 

e R (t) tte R , +5e R (t), (21) 

eiif) w e It0 + 5 ei (t), (22) 

where A = e^ j0 = e/ >0 = and A^o is one of the two roots of Eq. ffT6"|) . The linearized 
equations of motion are given by 

r -2A 

0-| 




" 8 A ' 




SN b 








5ei 





A -1 



£00-^ 





" 8 A ' 




8N b 




8e R 




8ej 



(23) 
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where the nonlinear interactions have been dropped and 

A gN-QN 



2 2v^Vo 

The eigenvalues of the Jacobian in Eq. (l23l) are given by 



2 

4 



1 ± y/l - 32eA/-/ 2 



and 



A 



3,1 



7 
4 



(24) 



(25) 



(26) 



1 ± VI - 8eA/ 7 2 j , 

respectively. Stability requires that the Jacobian has all eigenvaules with negative real parts, 
a condition which implies N > N/6, or equivalently, z 2 > 1/3. From this analysis we can 
conclude that for C < C CT i t , one root, z + > l/v^, is stable while the other root, z_ < l/\/3, 
is unstable. At C = C cr i t , these two roots merge and stability is completely lost as the 
system passes through a saddle node bifurcation. 




z , stable 



0.75 




0.25 



z , unstable 



0.1 0.2 0.3 0.4 
C 



FIG. 2: Stability plot for the regime with empty cavity, e = and \a = Xb = 0. The left panel 
shows Pz{z), whose roots give the stationary molecule fraction. The right pane shows the real 
roots of P${z), where the upper branch, z+, is stable for C < C cr u while the lower branch, z_, is 
unstable. 



As one can see from Eqs. (]14j) - (]15p , the factor g(N/2) corresponds to the maximal 
atom- molecule polarization that absorbs energy from the cavity mode (assuming iV 6 ~ N). 
In order to have e = stationary, this coherence term in Eq. f|T4l has to balance the driving 
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a) 




FIG. 3: Numerical solution of Eqs. ljl 21115]) . In a) the cavity field amplitude E(t) is plotted as 
a function of time for different values of the parameter C. the amplitude is scaled in units of 
&M = 2e/ , y. In b) the molecule fraction z 2 (t) is shown as function of time for different values of 
the parameter C. At C = C cr u ~ 0.38, the steady state disappears and large amplitude coherent 
oscillations between the atomic and molecular fields begin. This is accompanied by a nonzero 
cavity field. In these simulations Xa = Xb = W 6 in units of 7, small enough to be negligible in 
the regime C < Cent- 

term given by e. The constant C is a measure then of the ratio between these two terms, and 
when C > Cent, the driving is too strong to be compensated for by the absorption, and the 
steady state can no longer be maintained. Fig. ([3]) shows the cavity amplitude and molecule 
fraction as a function of time for different values of C. As one can see for C > C cr u, there 
are large amplitude nonlinear Rabi oscillations between the atomic and molecular fields. 
We focus next on the case of strong driving, i.e. C ^> C crit , which leads to a nonempty 
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cavity. Representing the cavity field using a polar representation, 

e(t) = E(t)e ie * {t \ (27) 
leads to the following dynamical equations 

M - QNh 

A = A Xa N - 2( Xb + 4Xa)N b + g b E sin(A + 9 e ), (28) 

2\J ^b 

N b = -gVN~b(N - 2N b )Ecos(A + 9 e ), (29) 
6 e = -| sin(6 e ) + ^^N~ b [N - 2N b ] sin(A + 6 e ), (30) 



and 



E = e cos(6 e ) 9 -\[W h \N - 2N b ] cos(A + 6 e ), (31) 



for the amplitudes and phases. Again we find the stationary values corresponding to A = 
N b = E = O e = in terms of the physical parameters. Here, that fact that E > and 
< N b < N/2 in the stationary state is used to assist in finding nontrivial solutions. First, 
for N b = 0, we get from Eq.( [29l that cos(A + B e ) = 0, which gives sin(A + O e ) = a = ±1. 
Using this in Eq.(|30l) gives, together with Eq.(fT71) 

sin(9 e ) = a| ^ z (1 - z 2 ) = a^z (l - , 2 ) . (32) 

In the regime of strong pumping C > 1 we have | sin(0 e )| < ^ <^ 1, which in turn implies 
cos(0 e ) ~ 1. These results when inserted in Eq. ([3lT) . yield the cavity amplitude 

E = — = e M , (33) 

7 

which is the maximal sustainable field amplitude and is what one would expect for a cavity 
with no atoms or molecules in the cavity. 

Using Eq. ( 1331) in Eq. ([281 , we finally arrive at a cubic polynomial for the molecule number 



z 3 - 3az 2 -bz + a = 0; (34) 

where 

a = aC- , (35) 

' (36) 



1 _ Xa. ' 

4 X6 
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The stationary molecule fractions thus only depend on the two parameters a and b for 
strong pumping. The solution manifold for Eq. (1341) . i.e. the surface of possible values of the 
molecule fraction < z(a, b) < 1 is shown in Fig. |H Here we see that there is a small region 
around < a < 0.5, < b < 1 where Eq. ([3~4"|) has two roots, which indicates the possible 
occurrence of bistability. 



FIG. 4: The real roots of Eq. [341 which gives the surface of stationary molecule fraction values as 
function of a and b for the regime with non-empty cavity, Eq > 0. Note that the left pane, (a), 
and the right pane, (b), show the same surface but for from different viewing perspectives. 

Next, we turn to the question of stability of the stationary solutions given by Eq.( l3"4"l) . 
Linearizing Eqs. fT251) - fl3Tl) for small perturbations around the stationary points 




a a 



e e (*)« e e , + 56 e (t) 



N b (t)nN + 6N b (t) 



A(t) w A + SA(t) 



E{t) ^E + 5E{t) 



(37) 
(38) 
(39) 
(40) 
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we get the Jacobian 



J 



B x B 3 

B 2 B 2 

n 7^3 _7 n 

u 4e 2 u 

f o o -I 



(41) 



where the matrix elements are given by 

B 2 
B 3 



- T fi(z), 



2eg N 

3 

"2./ ( f j 2 h(z) 



(42) 
(43) 
(44) 



V2\2 

and where we introduced the functions fi(z) = (1 + 3^ 2 )/2; 3 , /^(z) = z(l — z 2 ), and f$(z) = 
(1 — 3z 2 )/z, as well as the interaction constant x — 2(4% a — Xb)- The resulting secular 
equation for the eigenvalues is 



A 4 + A 3 + 



- + (e M -l)A-B 



A 2 + 



e M - 1 







r a 2 i 


B 


A - 


s + — 









0. 



where 



and 



5 



B2B3 



B\B 2 



(45) 



(46) 



(47) 



Using the Routh-Hurwitz stability criterion on Eq.( )45i) . we find that necessary conditions 
for the eigenvalues of the Jacobian to have non-positive real parts are 

1 



A > 



2(e M 



and 



and 



A 2 

eM 



> -B, 



-(e M -iy 

4 e M 



A 2 --(e M -l] 



k + ° 



A B < 

4 ~~ 



(48) 



(49) 



(50) 



The Routh-Hurwitz stability criterion is in general only necessary, and not always sufficient 



29| | . Here we find, however, that the two conditions fj4~9"]) and (1501) are sufficient as a numerical 
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check confirms that the roots of Eq. (1451) indeed have non-positive real parts when (1191) and 
fl50|) are fulfilled. In addition, we find that Eq(j39]) corresponds to loss of stability through a 
saddle-node bifurcation, whereas Eq. (f5T)|) corresponds to a Hopf-bifurcation as two roots of 
Eq. fl4"5l) cross into the right half of the complex plane. 
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FIG. 5: Region of stability in the A-B plane for different values of cm- The stability diagram 
was obtained from the Routh-Hurwitz criterion and explicitly checked by direct numerical solution 
of the roots of Eq. HSJ The region enclosed by the left and right borders represents the stable 
region where all roots have negative real parts. On the left boundary of this region there is a Hopf 
bifurcation and on the right boundary, a saddle node bifurcation. 



We first note that in general Eq. (l48j) is trivially fulfilled when the two inequalities (H9I) - 
( l50l) are fulfilled. Rewriting Eq.f l49l) in terms of the physical parameters we get 

For the case 7C 1, and a = — 1 combined with the assumption that the functions fa are 
limited, we have 

as an approximation for the stability boundary given by Eq. (l4"5j) . Fig. 5 shows the region 
of stability in the A-B plane. 

These results indicate that for specific choices of the parameters, both stationary states 
will be stable. As an example of the dynamics of the mean fields near the stationary states, 
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(52) 



200 




FIG. 6: Plot of time-dependence of the cavity field amplitude and molecule fraction for initial 
values, E = and z = 0.5. z = 0.5 was chosen to lie in the vicinity of the upper fold on the 
manifold of stationary solutions shown in Fig 4. E{t) is plotted in units of 7 and eu = 200. 

we solved Eq. (1341) for a « 3.3 x 10~ 4 , b=l/2 giving z ~ 0.5 and z ~ 10~ 7 . The values z = 0.5 
and z = were used as initial conditions in Eqs. (!28l) - (!3TT) . The results are shown in Figs. 
[6] and [7J for the cavity field and molecule fraction as functions of time. In Figs. [6] (a) and 
[7J (a), one sees that the cavity field amplitude approach cm exponentially for large times. 
The molecule fractions are shown in Figs. [6] (b) and [7J (b). For both cases, the solutions 
exhibit small amplitude oscillations around the stationary values, something characteristic 
of the presence of Hopf bifurcations. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the semiclassical dynamics of cavity assisted photo-association of 
molecules. Two different cases were described in some detail, namely weak and strong 
pumping. For weak pumping and zero two-body interactions, the steady state cavity field is 
zero. From Eqs. ffT4|) - ffl~5"]) we see that for a stationary state with e = 0, the pumping has to 
be exactly balanced by the atom-molecule coherence, which causes absorption. We interpret 
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FIG. 7: Same as Fig. 6 except that the initial value for the molecule fraction was chose to be 
z = 0, which lies close to the lower fold on the manifold of stationary solutions shown in Fig 4. 

this regime as that where the number of photons in the cavity is smaller than or equal 
to the maximum number of molecules. All of the photons injected into the cavity by the 
driving are used to convert a fraction of the atoms to molecules thereby leaving the cavity 
empty. When the pumping reaches a critical value, the system undergoes a saddle-node 
bifurcation and stability is lost. At this point the number of photons is greater than the 
maximum number of molecules and the system is able to undergo full amplitude nonlinear 
Rabi oscillations without fully depleting the cavity field. For strong pumping and non-zero 
two-body interactions, we find that there exists bistability in the molecule number in the 
prescence of a nonzero steady state cavity field. The bistable steady state molecule fractions 
are determined by the relative values of the atom-atom and molecule-molecule two-body 
interactions. The boundaries of the bistable region are marked by a saddle node bifurcation 
and a Hopf bifurcation. 

In a future publication we intend to go beyond the mean-field dynamics to explore the 
entanglement of cavity photons with molecular dimers. The form of H a tom-moi oc ihg(a^ 2 be — 
a 2 tfe*)/2 is of the same form as the four-wave mixing terms that give rise to optical phase 
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